Protected Qubits and Chern Simons theories in Josephson Junction Arrays. 
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We present general symmetry arguments that show the appearance of doubly denerate states 
protected from external perturbations in a wide class of Hamiltonians. We construct the simplest 
spin Hamiltonian belonging to this class and study its properties both analytically and numerically. 
We find that this model generally has a number of low energy modes which might destroy the 
protection in the thermodynamic limit. These modes are qualitatively different from the usual 
gapless excitations as their number scales as the linear size (instead of volume) of the system. 
We show that the Hamiltonians with this symmetry can be physically implemented in Josephson 
junction arrays and that in these arrays one can eliminate the low energy modes with a proper 
boundary condition. We argue that these arrays provide fault tolerant quantum bits. Further we 
show that the simplest spin model with this symmetry can be mapped to a very special Z2 Chern- 
Simons model on the square lattice. We argue that appearance of the low energy modes and the 
protected degeneracy is a natural property of lattice Chern-Simons theories. Finally, we discuss a 
general formalism for the construction of discrete Chern-Simons theories on a lattice. 



I. INTRODUCTION 

It is generally accepted that a quantum computer 
would have an enormous advantage over the classical one 
for the solution of many fundamental and practically im- 
portant problemsiS^ However, its practical implemen- 
tation presents a formidable challenge mostly because of 
the conflicting requirements posed by scalability and de- 
coupling from the environment. In particular, all scalable 
designs are based on the solid state devices but these 
are plagued by a strong decoherence. Quantitatively, in 
physics it is conventional to measure the decoherence by 
the quality ratio, Q, that is equal to the product of the 
decoherence time and a typical energy gap while in com- 
puter science one uses the error rate, R, defined as the 
probability of an error per time required for an individ- 
ual operation. In order to avoid excitations of higher 
energy states all operations should be performed slowly 
on the scale of the inverse energy gap, so R 3> 1/Q- 
The problem posed by the omnipresence of the decoher- 
ence in solid state devices is exacerbated by the fact that 
the error correction codes^j& require a small error rate 
for individual qubits (at the very least 10 -3 per logical 
operation that translates into the quality factors larger 
than 10 4 ) and lead to a huge increase in the number of 
qubits, effectively replacing one qubit by a lattice with 
L x L qubits with large l£>. Further, the efficient error 
correction requires that operations are done simultane- 
ously on all L x L qubits encoding a single error free bit 
which makes this scheme rather unrealistic. An alterna- 
tive would be to use individual qubits with a very high 
quality factor, much greater than 10 8 which would allow 
to perform calculations without (or with very little) er- 
ror correction. If the noise couples linearly to the energy 
difference between two states of the qubit representing 
and 1, this would mean that the physical noise should 



be less than 10~ 8 times than all other energy scales of 
the device. It is difficult to imagine a solid state physical 
system that is so well screened from the outside noise, in 
particular, it is difficult to imagine a Josephson junction, 
a Cooper box or SET where the motion of stray charges 
do not result in a significant 1/ / noise in electric or mag- 
netic fields or in the strength of the Josephson couplings. 
Further, a significant energy difference between the two 
states of a qubit results in a phonon emission^ that limits 
the quality factor of a typical Josephson device by 10 4 . 
So, the only logical possibility is that the physical noises 
do not affect the energy difference between lowest energy 
states in the lowest orders (one or more) in the strength 
of the noises. This is in principle possible because the 
effect of the physical noise is always represented by the 
sum of local physical operators (charge, current, etc). 
In the limiting case when the noise does not affect the 
energy difference between the two states in any finite or- 
der, these states form a protected subspace of the Hilbert 
spaced. Of course, the formation of such truly protected 
space becomes possible only for an infinite system. The 
implementation of these models in solid state (Josephson 
junction) devices was suggested i n 9 i 10 i n . From a prac- 
tical view point it is important to consider simpler but 
smaller systems which are protected from the noise in the 
given order, n, i.e. which are not affected by all physi- 
cal noises in all orders less than n. Clearly, any device in 
which two levels representing and 1 have a finite energy 
difference is susceptible to the fluctuations in the phys- 
ical quantity that sets this energy scale. For instance, 
in a Cooper box this would be the Josephson, Ej, and 
a charging energy, Eq of the individual Josephson junc- 
tion. Thus, even this limited protection can occur only 
when the two "working" levels are degenerate. 

It is well known that the stable degeneracy of the quan- 
tum levels is almost always due to a high degree of the 
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symmetry of the system. Examples are numerous: time 
inversion invariance ensures the degeneracy of the states 
with half-integer spin, rotational symmetry results in a 
degeneracy of the states with non-zero momentum, etc. 
In order for the degeneracy to be stable with respect to 
the local noise, one needs that the sufficient symmetry 
remains even if a part of the system is excluded. The sim- 
plest example is provided by the 6 Josephson junctions 
connecting 4 superconducting islands (so that each island 
is connected with every other) H In this miniarray all is- 
lands are equivalent, so it is symmetric under all trans- 
formations of the permutation group S±. This group has 
a two-dimensional representation and thus, pairs of the 
exactly degenerate states. With the appropriate choice of 
the parameters one can make these doublets the ground 
state of the system. The noise acting on one supercon- 
ducting island reduces the symmetry to the permutation 
group of the three elements which still has two dimen- 
sional representations. So, this system is protected from 
the noise in the first order (n = 2). The goal of this paper 
is to discuss designs giving the systems that are protected 
from the noise in the higher orders. Note that systems 
with higher symmetry groups, such as 5 junctions con- 
nected by 10 junctions (group S5) typically do not have 
two dimensional representations, so in these systems one 
can typically get much higher degeneracy but not higher 
protection. 

Generally, one gets degenerate states if there are two 
symmetry operations, described by the operators P and 
Q that commute with the Hamiltonian but do not com- 
mute with each other. If [P, ^ for any all 
states are at least doubly degenerate. Local noise term 
is equivalent to adding other terms in the Hamiltonian 
which might not commute with these operators thereby 
lifting the degeneracy. Clearly, in order to preserve the 
degeneracy one needs to have two sets (of n elements 
each), {Pi} and {Qi} of non-commuting operators, so 
that any given local noise field does not affect some of 
them; further, preferably, any given local noise should 
affect at most one Pi and Qi. In this case, the effect of 
the noise appears when n noise fields act simultaneously, 
i.e. in the n th order in the noise strength. Another im- 
portant restriction comes from the condition that these 
symmetry operators should not result in a higher degen- 
eracy of the ideal system. For two operators, P and 
Q that implies that [P 2 ,Q] = and [P, Q 2 } = 0. In- 
deed, one can construct the degenerate eigenstates of the 
Hamiltonian starting with the eigenstate, |0), of H and 
Q and acting on this state with P. The resulting state, 
|1) should be different from the original one because P 
and Q do not commute: [P, Q]^ 7^ for any ^. In 
a doubly degenerate system, acting again on this state 
with the operator P one should get back the state |1), so 
[P 2 ,Q] — 0. For a set of operators, the same argument 
implies that in order to get a double degeneracy (and not 
more) one needs that [PiPj,Q] — and [P, QiQj] = 
for any Indeed, in this case one can diagonalize 
simultaneously the set of operators {Qi}, {QiQj} and 



{PiPj}. Consider a ground state, |0), of the Hamilto- 
nian which is also an eigenstate of all these operators. 
Acting on it with, say, Pi we get a new state, |1), but 
since |1) cx (PPi)Pi|0) = P(PiPi)|0) cx P|0) all other 
operators of the same set would not produce a new state. 

In the rest of the paper we introduce models that pos- 
sess the symmetries satisfying these conditions (Section 
II), discuss their possible implementations in Josephson 
junction networks (Section III) and show that they are 
equivalent to the Chern Simons IL2 gauge theory (Sec- 
tions IV and V). Section VI summarizes our results. 

II. SPIN MODEL 

The conditions discussed at the end of the previous 
section are fully satisfied by the spin 5 = 1/2 model on 
a square n x n array described by the Hamiltonian 

Here a are Pauli matrices, note that the first term couples 
spins in same row of the array while the second couples 
them along the columns. It is not important for the fol- 
lowing discussion whether the boundary conditions are 
periodic or free, but since the latter are much easier to 
implement in a hardware we shall assume them in the 
following. Further, the signs of the couplings are irrele- 
vant because for a square lattice one can always change it 
by choosing a different spin basis on one sublattice. For 
the sake of argument, we assumed that the signs of the 
couplings are ferromagnetic, this is also a natural sign for 
Josephson junction implementations in Section III. The 
Hamiltonian JQ) was first introduced ini^ as a model for 
the anisotropic exchange interaction in transition metal 
compounds but its properties remain largely unclear. 

The Hamiltonian has two sets of the integrals of 
motion, {Pi} and {Qi} with n operators each: 

p* = IK 

i 

i.e. each Pi is the row product of of ■ while Qj is the 
column product of ct?-. Consider Pi operator first. It 
obviously commutes with the second term in the Hamil- 
tonian and because the first term contains two of ■ oper- 
ators in the same row, Pi either contains none of them or 
both of them and since different Pauli matrices anticom- 
mute, Pi commutes with each term in the Hamiltonian 
(Q. Similarly, [Qi, H] — 0. Clearly, different Pi commute 
with themselves, P 2 = 1 and similarly [Qi,Qj] = and 
Q 2 = 1, but they do not commute with each other 

{Pi,Qj} = (2) 
[PuQj} 2 = 4 



3 



so [Pi,Qj]\$>) ^ for any |\fr), thus in this model all 
states are at least doubly degenerate. Further, because 
PiPj contains two of ■ in any column, such product com- 
mutes with all Qk operators and similarly [QkQi, Pi] = 0. 
Thus, we conclude that in this model all states are dou- 
bly degenerate, there is no symmetry reason for larger 
degeneracy and that this degeneracy should be affected 
by the noise only in the n th order of the perturbation 
theory. 

To estimate the effect of the noise (which appears in 
this high order) one needs to know the energy spectrum of 
the model and what are its low energy states. All states 
of the system can be characterized by the set, {A^ = ±1} 
of the eigenvalues of Pi operators (or alternatively by the 
eigenvalues of the Qj operators). The degenerate pairs of 
states are formed by two sets, {Xi} and {— A^} and each 
operator Qj interchanges these pairs: Qi{\i} = {—Xi}. 
We believe that different choices of {Xi = ±1} exhaust 
all low energy states in this model, i.e. that there are ex- 
actly 2™ low energy states. Note that this is a somewhat 
unusual situation, normally one expects n 2 modes in a 
2D system and thus 2™ low energy states. The number 
2™ low energy states is natural for a one dimensional sys- 
tem and would also appear in two dimensional systems if 
these states are associated with the edge. Here, however, 
we can not associate them with the edge states because 
they do not disappear for the periodic boundary condi- 
tions. We can not prove our conjecture in a general case 
but we can see that it is true when one coupling is much 
larger than the other and we have verified it numerically 
for the couplings of the same order of magnitude. We 
start with the analytic treatment of the J z 3> J x case. 

When one coupling is much larger than the others it 
is convenient to start with the system where these others 
are absent and then treat them as small perturbations; 
in the limit J x = all columns are independent and 
the ground state of each column is fsing ferromagnet. 
The ground state of each column is doubly degenerate: 

i 1 )^ - ni t>« and i2>j = ni i>« s ivin s us 2 " d °- 

i i 

generate states in this limit. Excitations in each column 
are static kinks against the background of these states, 
each kink has energy 2J Z . Including now J x coupling, we 
see that it creates two kinks in each of the neighboring 
columns increasing thereby the energy of the system by 
8 J z so the lowest order of the perturbation theory is small 
in J X /8J Z . The splitting between the 2™ states occurs 
due to the high order processes which flip all spins in two 
columns. In the leading approximation one can calcu- 
late the amplitude of this process ignoring other columns. 
Thus, for this calculation we can consider the model with 
only two columns that can be mapped onto a single Ising 
chain in the transverse field in the following manner. The 
ground state of two independent columns belongs to the 
sector of the Hilbert space characterized by all PiPj = 1, 
it is separated from the rest of the spectrum by the gap of 
the order of 2J Z . Further, the Hamiltonian does not mix 
this sectors with different Pj, so in order to find the low 



energy states, it is sufficient to diagonalize the problem 
in the sector Pi = 1 . In this sector only two states are 
allowed in each row: | ||) and | J. J,), in the basis of these 
states the Hamiltonian is reduced to 

tfcol = -2J Z ]T T* t/ +1 -J x ^ (3) 

i i 

where r are Pauli matrices acting in the space of | ft) 
and | 11) states. This leads to the splitting 2A w 
(J x /2J z ) n (2J z ) (see Appendix A for the details of this 
calculation) between the symmetric and antisymmetric 
combinations of the two ferromagnetic chains in this 
problem. Thus, we conclude that the effective Hamil- 
tonian of the low energy states in the full system is 

3 

where r are Pauli matrices acting in the space of |1) and 
1 2) states describing the global state of the whole column. 
This effective low energy model also describes a ferro- 
magnetic chain in which the excitations (static kinks) 
are separated from the degenerate ground state by the 
gap 2A. In the basis of these 2" low energy states the 
operators Qj = tJ . 

We conclude that in the limit J z » J x 2" low energy 
states form a narrow (of the order of A) band inside a 
much larger gap, J Zl characterized by different eigenval- 
ues of Qj operators and by one value PiPj = 1. In the 
opposite limit J x 3> J Z low energy states form a narrow 
band inside the gap of size J x characterized by differ- 
ent eigenvalues of Pj operators and by the same value 
QiQj = 1- Consider now the effect of a weak noise in the 
former limit. To be more specific we consider the effect 
of the additional single site fields: 

H n = X! ^t^tj + hfjcrlj 

i,3 

The first term shifts (up or down) the energies of each 
ferromagnetic column by H* = /if j while the second 
term gives the transitions between up and down states 
in each column. These transitions appear only in the 
order n of the perturbation theory, so their amplitude is 
exponentially small: iJJ = (Y\_hfj/J z )J z . Thus, when 

i 

projected onto the low energy subspace this noise part 
becomes 

3 

The effect of the first term on the ground state degener- 
acy appears in the n th order of the perturbation theory 
in Hj / A and so it is much bigger than that of the second 
term because A becomes exponentially small asn^oo 
for J x <C J z . Note that although the effect of the /if -of j 
noise appears only in the large order of the pertubation 
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theory, it is not small because of the small energy de- 
nominator in this parameter range. Similarly, we expect 
that in the opposite limit, J z <C J x , the low lying states 
are characterized by the set of eigenvalues of Pi oper- 
ators, the effect of the hf jcrf j grows rapidly while the 
effect of the hfjafj noise decreases with the J x increase. 
We conclude that in the limits when one coupling is much 
larger than another ( J z 3> J x or J z <C J x ), the gap closes 
very quickly (exponentially) and the non-linear effect of 
the appropriate noise grows rapidly with the system size. 
These qualitative conclusions should remain valid for all 
couplings except a special isotropic (J x = J y ) point un- 
less the system undergoes a phase transition near this 
point (at some J x j J y = j c ~ 1). 
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FIG. 1: Energy spectrum of the 5x5 and 4x4 systems in the 
units of J z coupling as a function of J x / J z . We show energies 
of the lowest 40 states for 5x5 (upper pane) and lowest 20 
states for the 4x4 system (lower pane) . One clearly sees that 
at large anisotropy a well defined low energy band is formed 
which contains 2 5 states for 5x5 system and 2 4 states for 
4x4 one. In order to verify that low energy states are in one- 
to-one correspondence with P; eigenvalues for large J x /J z we 
have calculated the second lowest eigenstate in Pi = 1 sector 
(first one is the ground state). As shown in the lower pane, 
this state indeed has a large gap for J x /Jz > 1.2 

In order to check these conclusions we have numerically 
diagonalized small spin systems containing up to 5 by 5 
spins subjected to a small random field 7i* • flatly dis- 
tributed in the interval (—8/2,5/2). We see that indeed 
the gap closes rather fast away from the special J x = J y 
point (Fig. ^) but remains significant near J z = J x point 



where it clearly has a much weaker size dependence. In- 
terestingly, the gap between the lowest 2" states and the 
rest of the spectrum expected in the limits J z J x 
or J z <C J x appears only at J x /J z > j c with a practi- 
cally size independent j c ~ 1.2. We also see that the 
condition Pj = 1 eliminates all low lying states in the 
Jz <C J x limit where the lowest excited state in p = I 
sector is separated from the ground state by a large gap 
and in fact provides a lower bound for all high energy 
states. The special nature of this state appears only at 
Jx/Jz > jc- Clearly, the system behaves quite differently 
near the isotropic point and away from it but the size lim- 
itations do not allow us to conclude whether these differ- 
ent regimes correspond to two different phases (with the 
" isotropic" phase restricted to a small range of parame- 
ters j~ l < J x / J z < jc) or it is a signature of the critical 
region which becomes narrower as the size increases. Al- 
though we do not see any appreciable change in j c with 
the system size, our numerical data do not allow us to 
exclude the possibility that j c tends to unity in the ther- 
modynamic limit. We conclude that numerical data favor 
intermediate phase scenario. In contrast to this, the an- 
alytical resuls for two and three chains indicate that the 
transition occurs only at J z = J x point. Namely, both 
two and three chain models with periodic boundary con- 
ditions in the transverse direction can be mapped onto 
solvable models with transition at J z = J x : in the case of 
two chains the problem is mapped onto the exactly solv- 
able Ising model in transverse field as described above 
while the three chain model is mapped (see Appendix B) 
onto the 4 states Potts model in a similar way. The lat- 
ter is not exactly solvable in the whole parameter range 
but it obeys the exact duality that allows one to deter- 
mine its critical pointjMii£ , further, its exponents can 
be determined from the conformal field theory^ 6 . The 
mapping of three chain problem onto the Potts model is 
possible because the number of states of three spin rung 
for a given value of the conserved P operator is 4 while 
the number of different terms in the Hamiltonian that 
couples the adjacent rungs is 3. For a larger number of 
chains the number of states in each rung grows expo- 
nentially while the number of terms in the Hamiltonian 
grows only linearly making such mappings impossible. In 
this sense two and three chain models are exceptional and 
it is fairly possible that the intermediate phase appears 
in models with larger number of chains. 

Finally, we checked the effect of the hfjaf j on the 
ground state degeneracy splitting, our results are shown 
in Fig. |2 We see that, as expected, this disorder becomes 
relevant for J x <C J z while in the opposite limit its effect 
quickly becomes unobservable. We conclude that at (and 
perhaps near) isotropic point, the gap closes slowly with 
the system size and the effect of even significant disorder 
(5 = 0.1) becomes extremely small for the medium sized 
systems. 

Although it is not clear how fast the gap closes in ther- 
modynamic limit (if it closes exponentially fast the sys- 
tem never becomes truly protected from the noise be- 
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FIG. 2: Ground state splitting by random field in ^-direction 
for 5x5 and 4x4 systems. The random field acted 
on each spin and was randomly distributed in the interval 
(—0.05,0.05). Note that the effect of the random field in z- 
direction becomes larger for J x <C Jz as expected (see text). 
Because near J x = J z isotropic point the gap for 5x5 sys- 
tem is significantly smaller than the gap for the 4x4 system, 
this relatively large disorder has almost the same effect on 
these systems at J x ~ J z . We have verified numerically that 
decrease of the disorder by a factor of two leads to a dramat- 
ically smaller effects for 5x5 system confirming the scaling 
Ei — Eo oc 5 n discussed in the text; for J x / J z > 1.2 the 
difference Ei — Eo is difficult to resolve numerically. 

cause the effect of the high order terms might get very 
large), our numerical results clearly indicate that medium 
size (4 by 4 or 5 by 5) systems provide an extremely good 
protection from the noise suppressing its effect by many 
orders of magnitude. This should be enough for all prac- 
tical purposes. Further, if it is possible to construct the 
systems where PiPj = 1 (in other words with an addi- 
tional term in the Hamiltonian Hp = — A PiPj with 
significant A), this would eliminate the dangerous low 
energy states, leading to a good protection for all cou- 
plings strengths J x ^ J z . Indeed, in this case, we can 
repeat the previous analysis and conclude that the ef- 
fects of the noise appear only in the n th order and that 
now the perturbation theory in the "dangerous" noise, 
H* , implies the expansion in iJJ/A where A is no longer 
exponentially small but is the coefficient in the Hamil- 
tonian Hp. Thus, in this case these higher order terms 
become small. 



III. JOSEPHSON JUNCTION 
IMPLEMENTATIONS 



The basic ingredient of any spin 1/2 implementation 
in Josephson junction array is the elementary block that 
has two (nearly degenerate) states. One of the simplest 
implementation is provided by a four Josephson junction 
loop (shown as rhombus in Fig. |3Ji) penetrated by mag- 
netic flux $ /2ii2ii^. Classically, this loop is frustrated 
and its ground state is degenerate: it corresponds to the 
phase differences ±7r/4 across each junction constituting 



the loop. Two states (spin "up" and "down") then cor- 
respond to the states with phase difference ±7r/2 across 
the rhombus. For a isolated rhombus a non-zero (but 
small) charging energy, Ec = e 2 /2C, would result in the 
transitions between these two states with the amplitude 

r a E^E^e-*^^ (4) 

thereby lifting this degeneracy. Here s is numerical coef- 
ficient of the order of unity that was found iniS, s ~ 1.61 
and Ej is the Josephson energy of each junction. 

A. Simplest Josephson junction array. 

We begin with the Josephson junction array that has 
two sets of the integrals of motion, {Pi} and {Qi} dis- 
cussed above, which is shown in Fig. |31 This array con- 
tains rhombi with junctions characterized by Josephson 
and charging energies Ej > Ec and a weaker vertical 
junctions characterized by the energies ej, ec- As we ex- 
plain below, although this array preserves the integrals 
of motion, {Pi} and {Qi}, it maps onto a spin model 
that differs from l(T)l: we consider a more complicated 
arrays that are completely equivalent to spin model 
in the next subsection. The state of the system is fully 
characterized by the state of each rhombi (described by 
the effective spin 1/2) and by the small deviations of the 
continuous superconducting phase across each junction 
from its equilibrium (classical) values. Ignoring for the 
moment the continuous phase, we see that the potential 
energy of the array shown in Fig. |21is given by 

H z = -£J/ J Tj j T? +Ui , rf tj = JJ (5) 

i,j k<j 

Physically, the variable r' 3 - describes the phase of the 
rightmost corner of each rhombi with respect to the left 
(grounded) superconducting wire modulo ir. The right 
superconducting wire (that connects the rightmost cor- 
ners of the rhombi in the last column) ensures that the 
phase differences along all rows are equal. In the limit 
of a large phase stiffness this implies that the number 
of the rhombi with the phase difference 7r/2 should be 
equal for all rows modulus 2. This constraint does not 
allow an individual rhombus flips, instead a flip of one 
rhombus should be always followed by a flip of another 
in the same row. If, however, the phase stiffness is low, 
the flip of one rhombus can be also compensated by the 
continous phase deformations in the other rhombi con- 
stituting this row, we derive the conditions at which we 
can exclude these processes below. The simplest allowed 
process is the simultaneous flip of two rhombi in one row 

H x = -J2X(j-k)af tj af tk (6) 

where J x {k) is the amplitude to flip two rhombi a distance 
k apart. Both potential (J5J and kinetic © energies com- 
mute with the integrals of motion, {Pi} and {Qj}, so that 
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we expect that the main feature of this model, namely, 
the existence of the protected doublets will be preserved 
by this array. 




FIG. 3: Schematics of the array equivalent to the spin model 
with the interaction © in the vertical direction, (a) the main 
element of the array, the superconducting rhombus frustrated 
by magnetic flux <&o/2. Josephson energy of each rhombus is 
minimal for = and 8 = —n/2. Significant charging energy 
induces the transitions 8 = 0^8 = —n/2 between these 
energy minima, (b) The array geometry. The superconduct- 
ing boundary conditions chosen here ensure that PiPj = 1 
thereby eliminating all low lying states in the appropriate 
regime (see text), (c) The requirement that continuous phase 
does not fluctuate much while the discrete variables have large 
fluctuations is easier to satisfy in a very big arrays (L > 20) 
if one replaces the vertical links by the rhombi with junctions 
with ej, e~c frustrated by the flux 5>o(<5 + 1/2) 



As explained in the previous section, in order to 
achieve a really good protection one needs to eliminate 
all low energy states (except for the degenerate ground 
state) characterized by different values of the {Pi} and 
{Qj} operators. The array shown in Fig. |2]has a bound- 
ary conditions implying PiPj — 1 for any i,j because in 
this array the sum along each row of the phases across 
individual rhombus should be equal for all rows. Thus, 
for a sufficiently large tunneling amplitude J x (k) this ar- 
ray should have two degenerate ground states separated 
from the rest of the spectrum by a large gap. Physically, 
these two states correspond to two different values of the 



phase difference along each row. The quantitative con- 
dition ensuring that tunneling amplitude J x (k) is large 
enough depends on the range of J x {k). The simplest sit- 
uation is realized if only the nearest neighboring rhombi 
flip with the significant amplitude, J x . Because flip of 
the two nearest rhombi is equivalent to the flip of the 
phase on the island between them, in this case the spin 
model iFiKill is equivalent to the collection of independent 
vertical Ising chains with Hamiltonian 



H 



/ T X ■ 

x i>3 
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For J x ^> 2ej each chain described by this Hamiltonian 
has a unique ground state separated by the A = 2 J x from 
the rest of the spectrum. As the ratio ej/J x grows, the 
gap decreases. 

In the opposite limiting case of a very long range 
J x (k) = J x , one can treat the interaction © in the mean 
field approximation 



H x — JxL-r 



(7) 



At large J x the ground state of this system is also a dou- 
blet (characterized by (crf k ) = ±1) with all other excita- 
tions separated by the gap A = 2L X J X from the rest of 
the spectrum. As we increase the vertical coupling, ej, 
the gap for the excitations gets smaller. At very large 
ej the Hamiltonian is dominated by the ferromagnetic 
coupling in the vertical direction, so in this regime there 
are many low energy states corresponding to two possible 
magnetizations of each column. The magnitude of ej for 
which the gap decreases significantly can be estimated 
from the first order correction in ej. The dominant con- 
tribution comes from the transitions involving rhombi of 
the outmost rows. They occur with amplitude ej and 
lead to the states with energy A, so we expect that as 
long as ej < A, the system has a doubly degenerate 
ground state separated from the other states by gap of 
the order of A. 

The amplitude, J x (k), for the simultaneous flip of two 
rhombi can be found from the same calculation that was 
used iniS to calculate a single rhombus flip and the si- 
multaneous flip of three rhombi. If ec ^> Ec, the con- 
tribution of the vertical links to the total kinetic energy 
of the superconducting phase is small and can be treated 
as a small perturbation, in this case 



J x {k) a ETJ^E, 



3 / 4 pl/4 -2VBj/£c(l+cfc^) 



(8) 



where c ~ 1. Here the factor 2 in the exponential ap- 
pears because in this process one changes simultaneously 
the phases across two neighboring rhombi. Note that al- 
though the relative change in the action due to vertical 
links is always small, their contribution might suppress 
the flips of all rhombi except the nearest neighbors if 
EjEq/cq ^> 1. Note that even a relatively large ec (so 
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that Ec/ec "C 1) can be sufficient to suppress the pro- 
cesses involving non-nearest neighbors. We conclude that 
the low energy states become absent as long as 

ej < L eff J x (9) 

where L e // = 1/2 if EjEc/efj > 1 and L e ff w 
min(ec/v A ^j r Sc, if EjEc/e 2 c -C 1. These estimates 
assume that the main contribution to the capacitance 
comes from the junctions and ignores the contribution 
from the self-capaciatnce. If the self-capacitance is sig- 
nificant, the processes involving more than one island 
become quickly suppressed. 

We now consider the effect of the continuous fluctua- 
tions of the superconducting phase. Generally, a finite 
phase rigidity allows single rhombus flip, described by 
the ■ tafj term in the effective spin Hamiltonian. This 
term does not commute with the intergrals of motion 
Pi and thereby destroys the protected doublets. How- 
ever, for a significant phase rigidity the energy of a state 
formed by a single rhombus flip, U s f, is large. If, further, 
the amplitude t of these processes is small: t <C U s f, 
the states corresponding to single flips can be eliminated 
from the effective low energy theory and the protection 
is restored. If t > U s f, the protection is lost. 

We thus begin our analysis of the effects of the finite 
phase rigidity with the consideration of the dangerous 
single rhombus flips. Generally, the continuous phase 
can be represented as the sum of two parts: the one that 
it is due to the vortices and the spin-wave part which 
does not change the phase winding numbers. As usual 
in XY systems, it is the vortex part that is the most rel- 
evant for the physical properties. In particular, in these 
arrays it is the vortex part that controls the dynamics 
of the discrete subsystem. Notice that, unlike the con- 
ventional arrays, the arrays containing rhombi allow two 
types of vortices: half-vortices and full vortices because 
of the double periodicity of each rhombi. The flip of the 
individual rhombi is equivalent to the creation of the pair 
of half vortices. If the ground state of the system con- 
tains a liquid of half- vortices, these processes become real 
and the main feature of the Hamiltonian, namely the ex- 
istence of two sets of anticommuting variables is lost. We 
now estimate the potential energy of the half vortex and 
of the pair associated with single flip, U s f, and amplitude 
to create such pair, t. We begin with the potential en- 
ergy which is different in different limits. Let us consider 
the simpler limiting case when rhombi flips do not affect 
the rigidity in the vertical direction, it remains ej. Fur- 
ther, we have to distinguish the case of a very large size 
in horizontal direction and a moderate size because the 
contribution from the individual chains can be domini- 
nant in a moderate system if Ej ej. In a very large 
system of linear size L with rigidity ej in the vertical 
direction the potential energy of one vortex is 

E v = ny/Ejej ln(L) (10) 
while the energy of the vortex-antivortex pair at a large 



distance R from each other is 

U v {R)=^Ejej\n(R) (11) 

These formulas can be derived by noting that at large 
scales the superconducting phase changes slowly which 
allows one to use the continuous approximation for the 
energy density: E — ^Ej(d x 4>) 2 + \ej(d z (j)) 2 . Rescal- 
ing then the ^-coordinate by x — > xyEjJej we get an 
isotropic energy density E = ^^/Ejej(\7(j)) 2 . The con- 
tinuous approximation is valid if both rescaled coordi- 
nates x, z > 1. Thus, in a system with Ej ~ ej the 
formulas UlUlllfl remain approximately correct even at 
small distances R ~ 1, so a flip of a single rhombus 
creates a halfvortex - anti-halfvortex pair with energy 
E p rj Ej but the formulas become parametrically differ- 
ent in a strongly anisotropic system. Consider first the 
limit ej = 0. Here the chains of rhombi are completely 
decoupled and the energy of two half vortices separated 
by one rhombus in the vertical direction (the configura- 
tion created by a single flip) is due to the phase gradi- 
ents in only one chain, U^) = it 2 Ej/(2L), which appear 
because the ends of the chain have the phase fixed by 
the boundary superconductor. A very small coupling be- 
tween the chains adds J7 s y = (2/7r)e jL to this energy, 
so the total potential energy of the single flip inside the 
array is 

U sf = tt 2 Ej/{2L) + {A/n)ejL L « y/Ej/ej (12) 

This formula is correct as long at the second term is 
much smaller than the first one; they become comparable 
at L = \J Ej/ej and at larger L the potential energy 
associated with the single flip saturates at 

U sf = 1^/EjTj L » y/Ej/ej (13) 

where 7 « 3.3. Qualitatively, a single flip leads to the 
continuous phase configuration where phase gradients 
are significant in a narrow strip in x-direction of the 
length J Ej Jej and width ~ 1. The phase configu- 
ration resulting from such process is shown in Fig. 0] 
These formulas assume that the rigidity of the super- 
conducting phase remains ej which is, strictly speak- 
ing, only true if the discrete variables are perfectly or- 
dered in the vertical direction. Indeed, the coupling in 
the vertical direction contains ej cos(<^)r^ -rf +1 a which 
renormalizes to ej cos(<^)(r^ -t/ +1 ■) in a fluctuating sys- 
tem. In the opposite limit of strongly fluctuating rhombi, 
the average value of (T?jT? +1 j) becomes small, we can 
estimate it from the perturbation theory expansion in 
ej which sets the lowest energy scale of the problem: 
( T ij T i+ij) ~ e -j/(L e ffJ x ) which renormalizes the value 
of ej: 



This renormalized value of ej should be used in the es- 
timates of the vortex energy l|12f> or l|13fl ■ This does not 
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affect much the estimates unless the system is deep in 
the fluctuating regime. 



Two flips / 




A 

One flip \ 
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FIG. 4: Phase variation along the horizontal axis after a flip 
of a single rhombus (solid curve) and after a consequitive flips 
of two rhombi located at a distance twice the core size of each 
rhombi, Ej/ej. The horizontal axis shows the distance, 
X = y/ej/Ejx measured in the units of the vortex core size. 

Unlike potential energy, the single flip processes occur 
with the amplitude 

t = E^E^e-"^ 7 ^ (14) 

in all regimes. This formula neglects the contribution of 
the continuous phase to the action of the tunneling pro- 
cess. The reason is that both the potential energy HI3JI of 
the half- vortex and the kinetic energy required to change 
the continuous phase are much smaller than the corre- 
sponding energies of the individual rhombus, Ej, Ec- In 
order to estimate the kinetic energy, consider the contri- 
bution of the vertical links (horizontal links give equal 
contribution). There are roughly y/ Ej/ej such links, 
so their effective charging energy is about ec\/ej/Ej 
. If all junctions in this array are made with the same 
technology their Josephson energies and capacitances are 
proportional to their areas, so ej/Ej = Ec/ec = Vi m 
the following we shall refer to such junctions as simi- 
lar. In this case the array is characterized by two di- 
mensionless parameters, 77 1 and Ej/Ec S> I and the 
additional contribution to the charging energy, ^^E^ 1 , 
coming from vertical links is smaller than the one of the 
individual rhombi, Eq 1 and thus do not change the dy- 
namics. 

We conclude that the dangerous real single flip pro- 
cesses become forbidden if t <§; U s f where t is given by 
{T1J and U s f by {T^J) or {T3J). This condition is not dif- 
ficult to satisfy in a real array because amplitude t is 
typically much smaller than Ej. Further, for moderately 
sized arrays (with L = 5 — 10 which already provide a 
very good protection) the energy of a single rhombus flip 
is only slightly smaller than Ej, so the condition t <C U s f 
is not really restrictive. Note, however, that in order to 
eliminate low energy states of the discrete subsystem we 
also need to satisfy the condition © which implies that 



the tunneling processes should occur with a significant 
amplitude. While this might be difficult in the infinite 
array made from the similar junctions (with the same 
product of charging and Josephson energies), this is not 
really a restrictive condition for moderately sized arrays. 
One can choose, for instance, for a system of L x L rhombi 
with L = 5 — 10 Josepshon contacts with Ej = lQEc- 
This choice would give t ss 0.35Ec and J x « 0.2Ec for a 
system of disconnected horizontal chains. The condition 
t <C U s f is well satisfied. Choosing now the vertical links 
with e,j = Q.bEc and corresponding ec — 20Ec we get 
L e ff 1=3 5, so that the condition 10 is satisfied as well 
and there are no low energy states. It is more difficult, 
however, to eliminate the low energy states in the infinite 
array of coupled chains shown in Fig. and to satisfy 
the condition t <C U s f at the same time, especially if all 
junctions are to be "similar" in the sense defined above. 
This can be achieved, however, by replacing the vertical 
links by rhombi frustrated by the flux <f> = (1/2 + 5) <E>o 
with 8 <C 1 with each junction characterized by ej < Ej 
and ec ^ Ec- These rhombi would provide a signifi- 
cant rigidity to the continuous phase fluctuations (with 
effective rigidity e}) but only weak coupling (ej = 5ej) 
between discrete degrees of freedom. 

Finally, we discuss the effect of the finite phase rigid- 
ity on the amplitude of the two rhombi processes, J x {k). 
The condition that real single flip process do not occur 
does not exclude the virtual processes that flip conse- 
quitively two rhombi in the same chain. This would 
lead to an additional contribution to J x (k) ©. To es- 
timate this contribution we note that immediately after 
two flips the continuous phase has a configuration shown 
in Fig. 01 which is associated with the energy ~ U s f. 
The amplitude for such two consequitive flips is t 2 /U a f\ 
it can become of the order of J x (k) in a large system 
(where U s f is small). However, the amplitude of the 
full process involves additional action which further sup- 
presses this amplitude. This happens because the two 
consequitive flips lead to the high energy virtual state 
sketched in Fig. 0] and in order to get back to the low 
energy state the resulting continuous phase has to evolve 
dynamically. To estimate the action corresponding to 
this evolution, we note that its dynamics is controlled 
by \jEjjej junctions with charging energy ec- For the 
estimate we can replace these junctions by a single junc- 
tion with capacitive energy ec\/ e,j j Ej . Thus, the final 
stage of this process leads to the additional term in the 
action 5S ~ Ej/ec — r]{Ej / Ec)- Depending on the pa- 
rameter, \JEjEcjec = Tjy/ Ej j ' Ec, this additional con- 
ribution to the action is smaller or bigger than the total 
action, but even if it is smaller, it is still large compared 
with unity if Ej 3> ec- In this case, the processes that 
do not change the continuous phase dominate. We em- 
phasize again that in any case the transitions involving 
two flips in the same row commute with both integrals 
of motion P, Q and thus do not affect the qualitative 
conclusions of the previous section. For practically im- 
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portant similar junctions, it means the following. If 77 ^> 
yj Ec 1 1 Ej only nearest rhombi flip with the amplitude J x 
given by JHJ ■ If \/E c / Ej » 77 > E c / Ej the flips occur 
for the rhombi in the same row if they are closer than 
L eff < i 1 - 1 v / Ec/Ej. Finally for 77 < E c /Ej the dis- 
tance between flipped rhombi exceeds the size of the half- 
vortex and the two rhombi flips in a large (L 3> yEj/ej) 
array happen via virtual half-vortices in the continuous 
phase. 

In the discussion above we have implicitly assumed 
a superconducting boundary conditions such as shown 
in Fig. [3J These boundary conditions imply that in 
the absence of significant continuous phase fluctuations 
PiPj = 1. Physically, it means that if the array as a 
whole is a superconductor, it still has two states charac- 
terized by the phase difference, A<fi = or tt between the 
left and the right boundaries even in the regime where in- 
dividual phases in the middle fluctuate strongly between 
values and it. In this regime of strong discrete phase 
fluctuations, the external fields are not coupled to the 
global degree of freedom A(f> describing the array as a 
whole. In principle, it is also possible to have a similar 
array with open boundary conditions but in this case it is 
more difficult to eliminate low lying states because there 
is no reason for the constraint P t Pj = 1 in this case. 



B. Array equivalent to a spin model with local 
interactions. 

In order to construct the array equivalent to the spin 
model (pQ) we need to couple the rhombi in such a way 
that the transitions involving only one rhombus in a row 
are not allowed but the superconducting phase varies sig- 
nificantly between one rhombus and the next in the row. 
This is achieved if rhombi are connected in a chain by 
a weak Josephson link, characterized by Josephson en- 
ergy ej and Coulomb energy ec, so that Ej > ej and 
ec -C Ec, as shown in Fig. 03. In this case the simulta- 
neous tunneling of two rhombi which does not change the 
phase difference across the weak junction is not affected, 
its amplitude, J x , is still given by JBJl- 

In this array the flip of a single rhombus is due to two 
alternative mechanisms. The first one is that it involves 
the creation of the half-vortex (as discussed in the pre- 
vious subsection) and does not involve the change in the 
phase accross the weak junction. Alternatively, it can 
be due to the phase flip by n across the weak junction. 
Because ec <C Ec this process is slow and its amplitude 
is low: 

t = E z J A e^e- s '^l^ 

where s' » 1. Further, this process increases the en- 
ergy of the system by ej, so if its amplitude is small: 
t <C ej it can be completely neglected. Normally junc- 
tions with smaller Josephson energy, ej < Ej, have a 
bigger charging energy, not a smaller one as required 



here. To avoid this problem, we note that these weak 
junctions can be in practice implemented as a two junc- 
tion SQUID loops containing the flux $ = (1/2 + 5)$ 
as shown in Fig. The effective Josephson coupling 
provided by such loop is ej = sin(27r S)Ej (where E^p 
is the energy of the individual junction) while its effec- 
tive charging energy is ec — Ej/2. This allows to use 
bigger junctions for these weak junctions and has an- 
other advantage that it provides an additional control- 
ling parameter on the system. This construction seems 
somewhat similar to the partially frustrated rhombi that 
one needs to introduce as vertical links in large arrays 
discussed in previous subsection (Fig. 0;) but it serves 
a completely opposite purpose: it increases the density 
of the full vortices while keeping the discrete variables 
coupled. The partially frustrated rhombi, on the other 
hand, suppress the fluctuations of the continous phase 
while allowing independent fluctuations of the discrete 
variables in different rows. Such a dramatic difference is 
made possible by a combination of two reasons. First, 
the frustration has a different effect on these elements: 
in the case of the SQUID loops half flux eliminates the 
Josephson coupling completely while in a case of rhombi 
it leads to the exact degeneracy between two discrete 
states and to a rigidity for the continuous phase. Sec- 
ond, the values of the charging energies are rather differ- 
ent: in the case of SQUID loops the charging energy of 
its junctions completely dominates all processes in which 
phase changes across this loop thereby prohibiting the 
single rhombus flips. In the case of partially frustrated 
rhombi their charging energy is of the order of the charg- 
ing energy of the contacts in horizontal rhombi and thus 
it suppresses the double flips of distant rhombi in one 
chain but has relatively minor effect on the nearest ones. 

The J z interaction between the spins is provided by the 
pairs of the weak Josephson contacts as shown in Fig. 
with the Josephson energy ej < Ej and charging energy 
ec > Ec- These junctions do not affect the tunneling 
process of each rhombi but provide the weak interaction 
between them of the form |Q with the strength 

Jz = ej (15) 

Similar to the array discussed in the previous subsection 
we have to choose the parameters so that the energy of 
the half vortex and states resulting from a single flip is 
sufficiently high compared to the amplitude of the single 
flips. The discussion of the previous subsection carries 
over to this array. The only change is that the energy 
of the vortex in the infinite system contains the weakest 
link in the horizontal direction, i.e. ej instead of Ej: 

t < U sf 

U sf = vr 2 ej/(2i) + {A/^ejL L < y/ e.j/e., 
U s f = Jy/ejej L > y/ej/ej 

Finally, we have to ensure that the phases of the conse- 
quitive rhombi are decoupled and the interaction between 
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FIG. 5: Schematics of the array, (a) The elementary Joseph- 
son circuit emulating spin 1/2 consists of a four junction loop 
penetrated by magnetic flux $o/2. (b) Implementation of 
the spin chain by Josephson junction loops, here elementary 
rhombi are connected by weak links, the appropriate param- 
eters for these links can be obtained if each link in fact con- 
sists of two elementary Josephson junctions as shown in (c). 
(d). Full array implementing the spin model Q. The natural 
boundary conditions for this array are periodic, i.e. the point 
A should be connected to A', B to B' , etc. 



discrete degrees of freedom is purely local. This is sat- 
isfied if the continuous phase across weak junctions fluc- 
tuates strongly, i.e. that the energy of a usual (not half) 
vortex is smaller than the kinetic energy: y/ejej <C ec- 
Physically, it means that the array as a whole is an in- 
sulator, similar to the topological insulator considered 
iitli due to the full vortices that move in a vertical di- 
rection thereby decoupling different columns of rhombi. 
This condition does not contradict the condition t <C U s j 
because the latter involves the exponentially small am- 
plitude of flipping a single rhombus. If both conditions 
are satisfied the absolute value of the phase on each is- 
land constituting a rhombus fluctuates but the difference 
across the rhombus remains a slow varible taking two 
discrete values and it flips only simultaneously with the 
another one in the same row. Note that the interaction 
between discrete variables belonging to one column is 
due to the loop formed by these rhombi and two vertical 
junctions, it is therefore always local by construction, its 
value is given by ijTot . Repeating the arguments of the 
previous subsection, we see that in order to suppress the 
simultaneous flips of distant rhombi in the same row one 
needs also to satisfy the condition EjEc/^c ^ U but in 
contrast to the case of the regime discussed there, here cal terms, introducing an "electric" field (which is equal 



the conditions on the vertical junctions are not difficult to 
satisfy because one does not need to keep the long range 
order in a continuous phase. Under these conditions the 
dynamics of the array is described by the Hamiltonian 
©• 

Although in this regime the system as a whole is an 
insulator, it does not allow a half vortex to move across. 
So, physically, the two states of the global system can be 
observed in the array with the periodic boundary condi- 
tions shown in Fig[S] here two different states correspond 
to the half vortex trapped or not trapped inside the big 
loop formed by the array as a whole due to the periodic 
boundary conditions. 



IV. LATTICE CHERN-SIMONS GAUGE 
THEORIES WITH A FINITE ABELIAN GROUP 

In this section we discuss the general properties of the 
Chern-Simons theories with discrete gauge group Z n . 
Continuous Chern-Simons theories recieved a lot of at- 
tention in the past because they provide a natural mecha- 
nism for the flux attraction to the charged particles which 
thereby change their statistics^. In spite of some tech- 
nical difficulties, it is possible to construct lattice ver- 
sions of Chern-Simons theories22i2i, at least for continu- 
ous groups. The case of a discrete gauge group is slightly 
more delicate bcause one can not define canonically con- 
jugate pairs such as the "magnetic" (Aij) and the "elec- 
tric" ( . ) variables. Before we discuss the peculiarities 
associated with discrete groups, we briefly recapitulate 
the main properties of the continuous Chern-Simons the- 
ories with the non-compact and compact U(l) groups 
which we describe in the Hamiltonian formalism that we 
need to generalize these models to discrete groups. The 
U(l) Chern-Simons model is usually described by the La- 
grangian: 

£ = \KA 2 X + A*) - l -^B 2 + v{A x A y - A y A x ) (16) 

where B = d x A y ~ d y A x and dots stand for time- 
derivatives. We have used the gauge in which the time 
component A of the vector potential is zero. Because of 
this, we shall use only invariance under time- independent 
gauge transformations in this discussion. The canoni- 
cal variables conjugate to A x and A y are respectively 
IXj; = XA X + vA y and Tl y = XA y — vA x . The gauge 
transformations of the classical fields A p are the usual 
ones A p ~> A p + d p cf), but because of the Chern-Simons 
term this also induces a transformation of the conjugate 



fields II„ 



vt pa d a (i), where e xy 



1. 



In quantum theory, n p and A p become operators with 
the commutation relation [n p , Ap] = —i5 p ,p{r — r') and 
the gauge transformation is generated by the opera- 
tor R<p defined by i? = / d 2 rR(r)(f>(r) with R(r) = 
dpU p (r) + vB{r), since [R(r), A p (r')] = -id p S(r - r') 
and [R(r),Hp(r')] — —we pc7 d a 8{r — r'). In more physi- 



11 



to XA) in the classical theory) by E x — H x — vA y and 
E y = H y + vA x , the generator of the gauge transforma- 
tions can be expressed by R(r) = d p E p +2vB. It is simple 
to check that the gauge transformations commute among 
themselves and also commute with the Hamiltonian den- 
sity H — -rr\E 2 — ^/iB 2 . In the absence of matter, the 
equations of motion for the fields read: 



is the same as before reads: 



E n 



2v 

£ p a(-^-E„ 



[id a B) = 



XB = e pa d p E a 



(17) 



(18) 



This yields a massive branch of excitations, which are 
no longer purely transverse, but also develop a longi- 
tudinal component proportional to v. Their disper- 
sion relation is X 2 uj 2 — \v 2 + Xfik 2 . These propagat- 
ing modes are then pushed to very high energies in the 
limit where A <C v. In the presence of static exter- 
nal charges at positions r i: the Gauss constraint be- 
comes: R(r) — J2i SiS(r — rj). In the ground state, these 
charges induce a static field configuration according to: 
{2v - ^V 2 )5(r) = £i e ^(Y - n). So each particle is 
bound to a flux tube carrying a flux equal to f^- and 



smeared over a typical length £ 



(Am) 1 



The pure 



Chern-Simons limit is then recovered as £ goes to zero. 
In this limit, when one particle of charge e\ goes around 
an other charge , the total wave- function of the system 
is multiplied by an Aharonov-Bohm phase factor equal to 
exp(i^2-), so we a get a factor exp(i^f-) upon exchang- 
ing these two particles. We note that in the limit /i — > 
H commutes with operators Tl p = Tl p + ve m A, r These 
operators provide a generalization of the usual shift oper- 
ators (11^) for the Chern-Simons theory. Physically, the 
additional term in these operators appears because in the 
presence of Chern-Simons one can not only shift fields, 
one has also to change the phase of the wave function 
accordingly. 

We now turn to the lattice versions of the Chern- 
Simons theory. For the sake of simplicity, we shall work 
here on a square lattice, although these constructions 
could be generalized to other periodic systems. In this 
case, the vector potential describing the gauge field is de- 
fined on the links of the lattice, and will be denoted by 
Aij for the oriented link connecting sites i and j. If we 



reverse the orientation of the link, we obtain Aj 



-A, 



We shall adopt the Hamiltonian description from now on. 
In the absence of a Chern-Simons term, the local electric 
fields are simply the canonical conjugate operators Ely 
of the Aij J s. The corresponding commutation relations 
become: 

[A tj ,A kl ] = 

P«,n H ] = o 

[Aij, Ilfcz] = i&(ij),{M) 
The natural lattice Hamiltonian whose continuum limit 



H = ^ E n l + % + A * + 4w + A ^ 2 ( 19 ) 



<ij> 



ijkl 



In the presence of a Chern-Simons term, the electric field 
operators are modified. The electric field Eij will contain, 
besides Ely, terms associated to the vector potential in 
the direction perpendicular to the link (ij). On a square 
lattice, there are four links immediately perpendicular to 
this link, and containing either site i or j. Let us denote 
by Af(ij) this set of four links. To reflect the signs which 
appear in the continuous version discussed above, the 
links in Af(ij) have to be oriented in such a way as when 
(ij) runs from left to right, (kl) runs from bottom to top, 
as illustrated in Fig.[|Ja) below. 



Aj 



Aj 



(a) 



(b) 



FIG. 6: (a) An oriented link (ij) (dashed line) and the four 
oriented links adjacent to it (thick lines) which enter in the 
set J\f(ij). (b) For a given site i, the loop C(i) defined in text 
(thick lines). 

With these notations, the lattice Hamiltonian in the 
presence of a Chern-Simons term may be written as: 



H, 



cs 



<y> \ (kl)eM(ij) ) 



| Y^ A ^ + A J k + Ah + Ai ^ 

ijkl 



(20) 



The important fact is that the generalized electrical field 



operators Ei 



EL 



IE 



(ki)eM(ij) 



Aki no longer com- 



mute. The relation [E x (r),E y (r')] = —2iv8(r — r') from 
the continuous case becomes now: [Ey , Eki (r')] = — \v, 
whenever (kl) is one of the four links in Af(ij). Here (ij) 
is oriented along the positive x axis, and (kl) along the 
positive y axis. The apparent difference in the normal- 
ization of the delta function on the right hand sides of 
these expressions is compensated by the fact that M(ij) 
contains four elements. 

The energy spectrum for the eigenmodes associated to 
this system now reads: 

A 2 ^ = 4, 2 cos 2 (^)cos 2 (^)+4AMsin 2 (^)+sin 2 (^)) 
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This spectrum is specially interesting in the limit where 
(i goes to zero, since then it exhibits lines of zero modes 
at the boundaries of the first Brillouin zone, namely for 
k x = ±7r or k y = ±7r. These modes are directly related to 
two remarkable non-local conservation laws which appear 
in the limit of vanishing ji. More precisely, for each row 
and each column, we may define two conserved quantities 
P row and Qcoiumn in the following way: 

P row = J2 (-l) x ^ W +* (21) 

j £E row 

Qcolumn ^ ^ ( 1)^ P?,;/+y (22) 

j £ column 

Unlike the row and column operators discussed in Sec- 
tion II these operators mutually commute. Note that in 
the [i = limit, the system exhibits a large set of local 
conserved quantities Ily, i.e. one for each link, defined 
by: 

IT, = IT, + | Y*(ki)eM{ij) A ki- (23) 

Similar to the local electric fields, these variables 
do not mutually commute, so they cannot be si- 
multaneously diagonalized. Furthermore, they are 
not invariant under gauge transformations, since the 
electric field opeators are gauge invariant and since 
Eij -flij = -|Z)(fcOeJV(»j) A kl which is clearly gauge 
dependent. This situation is very similar to what occurs 
in the process of quantizing the motion of a particle on 
a plane in the presence of a uniform perpendicular mag- 
netic field. The local electrical field operators are analo- 
gous to the gauge invariant velocity operators, P — eA, 
in the particle problem. Because of the magnetic field, 
the two components of this vector no longer commute. 
Since the vector potential A is not translation invari- 
ant, the usual translation operators have to be combined 
with gauge transformations in order to commute with 
the kinetic energy (P — eA) 2 . These deformed genera- 
tors of translations are analogous to the two components 
of the vector II = P + eA. Apparently, the only way 
to construct gauge-invariant symmetry operators in the 
lattice gauge model is to use the non-local combinations 
P r0 w and Qcoiumn defined in (|21I22|I . Note that a model 
with similar conservation laws, for interacting Bosons on 
a square lattice, has been analyzed 

A next step is to construct a lattice Chern-Simons 
gauge theory for the continuous but compact U(l) group. 
This is simply achieved in the absence of a Chern-Simons 
term, by assuming that the vector potential variables Ay- 
are periodic, with a period chosen for instance equal to 
2tt. This implies that the spectrum of the conjugate op- 
erators Ilj,- are discrete, containing only integer values. 
Because of this, the form of the Hamiltonian has to be 
modified from equation (|19f) above, and a natural choice 
respecting the requirements of gauge invariance and pe- 
riodicity in the gauge potentials reads: 

H=^Y1 ^l-^^cosiA^+Ajk+Aki+Au) 2 (24) 

<ij> ijkl 



Using this Hamiltonian as a starting point we add a 
Chern-Simons term by the standard deformation of the 
electrical operators: 

Hcs = 5xH ( n ^4 £ Akl ) 

<ij> \ (fci)eAT(ij) / 

C0S ( A v + A Jk + Am + An) 2 (25) 

ijkl 

We now discuss what can be the Hilbert space associ- 
ated with this Hamiltonian. This is a non-trivial prob- 
lem because we can no longer impose the periodicity 
in the local gauge potentials in the usual way, assum- 
ing that the wave-function of the system considered as a 
function of the Ay's is 2ir periodic with respect to any 
Aij . This is not possible because the naive shift operator 
gnaive = e xp(-i27rIIy ) that changes Ay into Ay + 2tt 
no longer commutes with the kinetic part of the Hamilto- 
nian. The appropriate definition of these shift operators 
becomes Sy = exp(— i27rHj,), with the Hij defined in 
(|23|l . It order to construct invariant states under this full 
set of shift operators, we need them to be mutually com- 
muting. This is realized only for a discrete set of values of 
v = ^, where m is any integer. Thus, the compact gauge 
theory is compatible with the Chern Simons term only 
for special, "quantized", values of v (see alsoSiSi) This 
statement is very similar to the well-known fact that in 
order to quantize the problem of a particle on a torus in a 
perpendicular uniform magnetic field, the total magnetic 
flux through this torus should be an integer mutiple of 
the flux quantum - . In one particle problem this require- 
ment simply expresses the need for mutual commutation 
between two magnetic translations that are used to con- 
tract the torus from an infinite plane. 

Note that, even for these special values of v, the 
Hilbert space of the theory is peculiar because the shift 
operators Sy are not gauge invariant. More precisely, 
let us denote the generator of the gauge tranformation 
sending Ay into Ay + a, for any site j connected to 

site i by Ui(a) = cxp(— iocY^f' fty). In this expression, 
the sum is taken over the nearest neighbors j of site 
i. Clearly, these unitary operators commute with the 
Hamiltonian l|25|) . For a site i, and a link (jk) belong- 
ing to Li, where Hi is the oriented loop (see Fig. HJb)) 
defined by the edges of the square built from four ele- 
mentary plaquettes whose center is located at site i, we 
have the following relations: 

Ui(a)S jk = e ±im,a S jk U i {a) (26) 

The sign in the phase factor depends on the orientation 
of the link (jk), and is negative if it is oriented along 
Ci. This equation implies that one cannot enforce at the 
same time a condition of gauge invariance on the physical 
states, and invariance of the wave-function under the 2~k 
shift operators. If we choose to work in a basis of eigen- 
vectors for the gauge transformations, namely states \%) 
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satisfying Ui(a) = exp(iaQi)| v I'), where Qi is the electri- 
cal charge at site i, applying Sjk to for (jk) in Li 
modifies the charge Qi by ±m if f = — . Conversely, a 
subspace where Sjk\^) — e l9 i k Y$>) is not gauge-invariant, 
since after applying Ui(a) on such states, we get the new 
periodic conditions with 9[- k = djk + irva, for any link 
(jk) belonging to Li with the corresponding orientation. 
However, it appears that Ylk @jk = Sfc ®J k because 
rife ^ Sjk — Uj(2Tt) commutes with any gauge transfor- 
mation. This relation shows that the gauge invariant 
quantity Ylj? Ojk is nothing but 2wQj modulo 2n. 

The properties of this model have been investigated by 
several groups. In the absence of a Chern-Simons term, 
the periodicity of the U(l) gauge field allows quantum 
tunneling processes where the local flux on a plaquette 
changes by ±2ir. In a 2 + 1 dimensional path integral 
description, these instantons (called here monopoles) in- 
teract via a Coulomb-like 1/r potential, leading to De- 
bye screening in this monopole plasma. The proliferation 
of such tunneling events drives the system into a phase 
where the magnetic variables Aij are strongly disordered, 
and in which external static electric charges interact by 
a confining potential increasing linearly with their sepa- 
ration2SiS&. In the presence of a Chern-Simons term, the 
properties of the monopole plasma are deeply altered, 
and several groups have reached the conclusion that a 
linear interaction now binds pairs of monopoles of oppo- 
site charges, thus destroying the confinement of electrical 
chargea 21 i 27 i 28 . 

To extend this construction to a discrete Z„ group, we 
replace the continuous vector potential on the link join- 
ing sites i and j by = ^-Pij, where p^ is an integer. 
In the absence of a Chern-Simons term, the generator of 
the gauge transformation based at site i sending pjk into 

Pjk + Sji — Ski is Ui = JJ^ n+j, in analogy with the con- 
tinuous case discussed above. The unitary operator tt^ 
is analogous to the operator exp(— i—Hij) of the contin- 
uous model, namely it transforms into A^ + 2i:/n or 
equivalently p^ into pij + 1. In order to attach flux to 
particles, the generator Ui has to be modified by a phase 
factor exp(-i|(^) 2 X;( J fe) e £ i Pjfe), where the (jk)'s be- 
long to the loop Li already defined above (see Fig.^Jb)). 
For each value of is, we define the generators by: 

Ui=f[4^M-i^ f E »*) ( 27 ) 

3 Uk)dC(t) 

This definition preserves the fact that these generators 
mutually commute. Note that since the generators Ui 
commute simultaneously with the local fluxes, a conve- 
nient gauge-invariant basis of the Hilbert space is ob- 
tained by simultaneous diagonalization of the local fluxes 
and the f/j's. With a discrete basis to describe each link, 
the kinetic part in the Hamiltonian i|25|) has to be mod- 
ified. The most natural way to do this is to replace the 
local electrical field By by a gauge-invariant operator £ 



which shifts A^ by the minimal possible amount 2n/n. 
This operator is defined as follows: 

(ki)e^(ij) 

The Hamiltonian may now be written as: 

<*j> 

- ^E cos ( — (Py +Pjk +Pkl + Pli)f (28) 
ijkl 

where £~j is the Hermitian conjugate of To re- 

cover the Hamiltonian (|25|l in the limit where n becomes 
very large, we notice that acts in the same way as 
cxp(— i— Eij), therefore we have to choose A„ so that 
2(2 7 r/n) 2 A- 1 =A- 1 . 

With these notations, the operators shifting the link 
variables p^ by one unit, and which commute with all 
the gauge-invariant kinetic terms read: 

n± =7T±exp{-i~A 2 E Pw) 

since they are analogous to exp(— i-^-fljj) in the contin- 
uous model. The 2it shift operators Sij previously intro- 
duced are then equal to (7fy) n . Note that the parameter 
v is quantized, in the same way as before (namely v = ^ , 
with integer m) , since we impose the model to be periodic 
when p^ is changed into p^ +n. More precisely, as for the 
compact U(l) group, this notion of periodicity requires 
the mutual commutation between all the 2ir shift oper- 
ators S^. Models obtained from two values of v which 
differ by an integer multiple of 2n 2 /ir are clearly identical. 
We also note that changing m into 2n 2 — m amounts to 
replacing all the phase factors entering in the definition of 
various operators such as Ui by their complex conjugates, 
and the corresponding models exhibit similar properties. 
In the special case m = n 2 , the operators mutually 
commute, so this case is equivalent to m = 0. Therefore, 
it is sufficient to choose m in the set of integers between 
and n 2 — 1 . Among those n 2 possible values of v, there is 
an interesting subset of n values for which m is an integer 
multiple of n. If this condition holds, the generators of 
the elementary gauge transformations Ui defined in equa- 
tion H27fl commute with all the 2ir shift operators Sjk- It 
is then possible to apply simultaneously the condition of 
gauge invariance and the periodic boundary conditions 
Sjk\^) = 1^) for the allowed physical states. 

If v = 0, the resulting discrete gauge theory (without 
Chern Simons term) has two regimes according to the 
relative size of the two terms entering in Hcs- When 
fi ^> A^ 1 , H{y = 0) describes a phase where local- 
ized flux excitations, called fluxons, have the energy gap 
fi(l — cos( — )). In this phase, quantum fluctuations of 
the magnetic variables p^ are small, and the effective 
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interaction between two static external electric charges 
varies logarithmically with their separation. Notice that 
this phase owes its existence to the discrete nature of the 
symmetry group, as illustrated by the vanishing of the 
corresponding energy gap as n is taken to infinity . A 
small A^ 1 term simply gives some amount of dispersion 
to these excitations. When A is decreased further below 
a critical value A^, the fluxon gap closes, and the system 
enters the charge confining disordered phase similar to 
the one obtained in the compact U(l) case, for all values 
of A. 

Let us first consider the effect of a switching on a 
Chern-Simons term in the former regime, where the po- 
tential energy (proportional to /x) dominates. When A ra 
is very large, we do not expect the flux binding mecha- 
nism to operate. Indeed, a unit charge at site j corre- 
sponds to imposing the Gauss law constraint Uk\^) = 
exp(i^-Sj k )\^>). When A" 1 is small, it is energetically 
more favorable to keep a vanishing flux everywhere, be- 
cause of the low value of the kinetic term compared to 
the fluxon gap. So we expect the flux attachment mech- 
anism to work only if A„ is smaller than a critical value 
A*. When A„ is further reduced, below A£, the fluxon 
gap eventually closes, and the qualitative properties of 
the system change drastically. In Section V we present 
a simple perturbative estimate of the critical value AJj 
and argue that it is in fact equal to A*, i.e. both tran- 
sitions happen simultaneously. As already discussed for 
the case of the compact U(l) group, the presence of the 
Chern-Simons term modifies deeply the picture obtained 
for vanishing v in the strongly fluctuating regime of small 
A„. Note that, by contrast to the v = case, analysis 
of the A — > limit is difficult since the Eij operators no 
longer commute if they involve two perpendicular links 
sharing a common site. In the n — > oo limit, we expect 
to recover the continuous, but compact U(l) theory for 
which we still do not know how to write down explicitely 
the ground state wave- function, even in the A — » limit. 



V. MAPPING CHERN-SIMONS THEORIES 
ONTO SPIN MODELS 

Here we shall study in more detail the interesting limit 
of the vanishing magnetic energy and show explicitely 
how, in the n = 2 case, it maps precisely on the mod- 
els studied in the beginning of this paper. As we have 
discussed, it is possible to propose a design of Josephson- 
J unction arrays which directly implements this limit. 

As a first step, it is convenient to introduce a ba- 
sis in the Hilbert space of gauge invariant states (i.e. 
states such that Uil 1 ^) = for any site i), which 
keeps track of the flux variables. For any square pla- 
quette (ijkl) centered at r, the corresponding flux a r is 
the integer pij +Pjk + Pu +Pu- For any flux configura- 
tion {cr}, we may then define a gauge-invariant quan- 
tum state ^({ov})). Our main task is to represent the 
algebra of gauge invariant operators £^ in such a basis. 



These operators obey two families of constraints. First 
we have: 

= eX p(i2n^)£+£+ (29) 

Here the sign corresponds to the geometry when link (jk) 
is perpendicular to (ij) and located on its left. Notions of 
left and right are defined for an observer moving along the 
link in agreement with its orientation. We have used the 
condition v — m/ir. Second, these operators are related 
to the generators of local gauge transformations by: 

(i) W 

n £ io = cx p(* 27r ^ J2 a ^ ( 3 °) 

3 r 

As usual, the product in the left-hand side runs over the 
nearest neighbors of site i, whereas the sum in the right- 
hand side involves the four plaquettes adjacent to i. In 
this expression, attention should be paid to the ordering 
of the various operators. We assume that the two op- 
erators involving bonds along a given direction directly 
follow each other. Once this is enforced, any of the eight 
possible residual permutations compatible with this cri- 
terion does not change the result. For any oriented link 
(ij), let us call r (resp. r') the adjacent plaquette lo- 
cated at the left (resp. right) of (ij). The operator £^ 
decreases the local flux ay and increases o~ r by one unit. 
We see that £f- should be proportional to o-+a~, up to a 
phase factor which depends on the configuration of local 
fluxes on the whole lattice. This phase factor is required 
in order to satisfy the constraints (|29I3U[1 above. In a gen- 
eral case this phase factor might become a very non-local 
function of the flux configuration but it remains simple 
in the case of ~E-2 m = 2 model. We discuss now its 
construction in different cases, starting with the simplest 
ones. 

We start with the simplest case n — 2, and m = 2. 
Then, the electrical operators £^ on two adjacent and 
perpendicular links anticommute. Furthermore, accord- 
ing to equation (|2*H|) . the generators of the local gauge 
transformations commute with the shift operators Sjk = 
(tt^.) 2 - It is then possible to impose simultaneously the 
gauge invariance constraint and the ILi periodicity on 
the links. As a result, the two operators and <j~ are 
equal. We may then represent the above algebra by the 
following Pauli operators, associated to a quantum Ising 
model attached to the plaquettes of the lattice: 

£y = ofcy for vertical (ij) (31) 
or 

£y = o%o%, for horizontal (ij) (32) 

exp(i7rcr r ) = cr*/ (33) 

The symmetry operators defined in equations l|21|l and 
(|22|l are easily generalized to the case of a discrete group. 



15 



We now have: 

P =U £ ti+* ( 34 ) 

j^row(r) 

Q s = II Ef?;f (35) 

j£column(s) 

In these equations, the symbol sign(j) stands for + or 
— , according to the parity of Xj (resp. yj) in the first 
(resp. second) equation, row(r) denotes row number r 
and column(s) denotes column number s. For a generic 
(not boundary) row and column these operators com- 
mute. 

It is now possible to recover the symmetry operations 
Pi and Qj introduced in Section II for Chern-Simons the- 
ories with proper boundary conditions. First, notice that 
for 7Li model operators P r and Q s correspond to the prod- 
uct of two P (or two Q) operators: P r = P r P r +\ when 
expressed in terms of the spin operators of Section II. 
If operator Pq = 1 is trivial (Pq = 1), Pi = Pi and 
we can recover a single P r operator taking the product 
of all P r > operators with r' < r. In terms of the gauge 
theory, the Chern-Simons Hamiltonian (Eq. 1)28(1 1 in its 
/i = limit should not contain the dynamical variables 
(gauge fields) associated with the bonds along the edge 
of the lattice. In this case the operators Pi and Qi in- 
deed do not commute because they contain only one pair 
of non-commuting electric fields located in the corner of 
the lattice. With this assumption, we obtain two families 
of operators commuting with the Chern-Simons Hamil- 
tonian at fi = 0: 

p _ TT csign(j) 

Zj^column(s) 

Here sign(j) is defined as above. These operators satisfy 
the generalized commutation relations: 

fry 

P r Q s = exp(i27r— )Q a P T (36) 

In the special case m = 2, these operators anticommutc. 
Using equations l|31l) and (|32ll above, we may write: 

Pr= n < 

r£row(r) 

q s = n a r 

rEcolumn(s) 

We now discuss the mapping in a more general case. 
First, we assume that the generators of local gauge trans- 
formations still commute with the shift operators Sjk on 
the links. According to equation l|26ll . this requires m to 
be a multiple of n. In this case, we may still view the 
local fluxes a r as defined modulo n. This implies in par- 
ticular that (cr+) n can be chosen to act as the identity 



operator. After enforcing the gauge invariance and the 
2„ periodicity in the bond variables, we may represent 
S+ as: 

£+ = exp(a r , r '({a r })) ( T+ ( T r ; (37) 

The "statistical" gauge field A r .r' should not be con- 
fused with the original link variables Aij . This new 
entity is imposed to us by the necessity to satisfy the 
constraints i|29!3UI) . Note that as usual, there is a large 
amount of freedom in the definition of -4 r ,r', reflecting 
the arbitrariness in choosing a global phase for each state 
^({crr})). If we could ignore the first constraint 12911 . 
our system of fluxons would be completely equivalent to 
a collection of particles obeying fractional statistics, since 
the second constraint l|3UI) relates the total statistical flux 
seen by a fluxon hopping around an elementary plaquette 
of the dual lattice centered at i, to the number of fluxons 
in the immediate neighbor hoood of i. The presence of the 
first constraint is an original feature of fluxon dynamics 
in Chern-Simons theories. 

Let us now show how to construct explicitely one real- 
ization for this statistical gauge field A r , r ' ■ This amounts 
to making a definite choice for the global phase of the 
basis states ^({cr})). If we have one fluxon on the pla- 
quette centered at r, it is possible to represent this by 
a string-like configuration of bond variables p^ (r) such 
that vanishes on any horizontal bond, and on most 
vertical bonds, with the exception of all the links located 
on the same row as r and on its right, for which it takes 
the value 1. This may be summarized by the following 
definition: 

|r)=Pinv 1] 41°) 

(ij)>r 

Here, the operator Pi nv is the projector on the subspace 
of the gauge invariant states. The notation (if) > r 
stands for all the links (ij) on the right-hand side of the 
plaquette center r and on the same row, and the reference 
state |0) is simply the state where all the bond variables 
Pij are equal to 0. For a general flux configuration {er r }, 
we simply choose a reference configuration Pij({cr r }) ob- 
tained by superposing the configurations associated to 
each fluxon excitation in the system. Therefore, the in- 
tegers Pij({<J r }) are defined as follows. They are equal to 
zero for any horizontal bond (ij). For a vertical bond, 
we set: 

r<(ij) 

Using this, we define: 

|*(K}))=P inv n(^r 3({CTr}) |0) (38) 

(ij) 

We could have used operators 7r^ instead of £7r to 
generate a state with the same desired flux configura- 
tion {cr r }, but the advantage of £ f- is that they commute 
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FIG. 7: The electrical field operator, £i a i 1 moves one flux 
quantum from the plaquette r to r' . In order to represent 
this process in the gauge where all horizontal links Aij are 
0, we use the gauge transformation on the string of sites 
ii . . .i„ that produce the electric field operators on vertical 
links (shown by light arrow) and the phase factors associated 
with the shaded area. In the flux representation the effect of 
vertical electric field is to remove the flux from plaquette r and 
to add one flux to plaquette r' . Note that the phase factors 
associated with all plaquettes except r and r' axe doubled. 



with the projector Vi nv . These states |\&({(7 r })) form 
a complete basis of the gauge invariant subspace, since 
any configuration {pij} of bond variables producing the 
flux pattern {cr r } may be deduced from the configuration 
{PijiWr})} by a gauge transformation. With this choice 
of gauge, we have: 



E ij = CT r °"r 



(39) 



for any vertical link (ij). This means that A r y vanishes 
whenever the vector joining r and r' is equal to ±x. For 
a horizontal link (ij) oriented from left to right, we make 
a repeated use of the second constraint (|3l7|) above for all 
the sites k located on the same row as (ij) and on its 
right, which, if applied on a gauge invariant state, yields: 



= cxp ( i2n^(o- r + a r > +2 ^ <?>") ] Mi 



r" >r,r' 

where the "string" operator M.^ is defined by: 

m ^ = n s m n s ki 

(kl)>r (kl)>r' 

As before, r (resp. r') denotes the adjacent plaquette 
located above (resp. below) the oriented link (ij). The 
notation (kl) > r refers to all the vertical links (kl) on 
the right of plaquette r and on the same row. In our 
gauge the string operator becomes 

M^\{a r })=a+a-\{a r }) 

Finally, we get 



£ij = 0~r a r> CX P ( * 27r ^2'( ' r + °V + 2 ^ <7 r «) 



v" >r,r' 



where the formula applies only on the gauge invariant 
subspace. This shows that 



exp(iA r y) 



TYl 

exp I i2n—(a r 



r" >r,r' 



whenever the vector joining r and r' is equal to ±y. The 
right-hand side of this expression is most of the time 
highly non-local. The only exception is the case when 
2m is an integer times n 2 . For m chosen in the inter- 
val [0,n 2 — 1], this occurs only when 2m = n 2 , which 
requires n to be even. If this condition is satisfied, we 
get: 

exp(iA r , T +y) = exp(i7r(er r + oy)) 

In physical terms, this set of models corresponds to an 
effective Bosc statistics for fluxons. 

Let us now address the most general case. Since ac- 
cording to i|26|) . exchanging Ui and Sjk involves a phase 
factor exp(±«27r^), we may always diagonalize simulta- 
neously the local gauge generators E/j with the operators 
S™ k where n' is equal to n divided by the greatest com- 
mon divisor of m and n. This implies that dimension of 
the Hilbert space associated with one flux is increased 
to N — nn 1 and we may view the local fluxes as ele- 
ments of Zjy. All the expressions already written are 
then valid. In particular, we may check that exp(i„4 r . r /) 
is unchanged, if we change <r r into o~ r + Nb r , where b r is 
any integer. 

The representation in terms of fluxes constructed in 
this Section becomes especially convenient when the gap 
of a single fluxon is the largest energy scale in the prob- 
lem which occurs if /i >> A,^ 1 (see Hamiltonian (|28|l ). In 
this limit, the term proportional to A^ 1 induces tunnel- 
ing processes where a fluxon jumps from a plaquette to 
one of its neighbors. Using the constraints (|29fl l|3U|) . we 
see that the operator which moves a fluxon around the 
elementary plaquette of the dual lattice centered at site i 
is simply equal to Ui, provided that it acts on a state with 
only one fluxon on a plaquette adjacent to i. Indeed, this 
operator may be written as: ^-y E t^ E t+y E t+^ n0 " 
tice that this ordering is different from the equation l|3U|l . 
To recast this product in the form of the latter expres- 
sion, one may permute the last two operators on the 
right, which according to l|29() produces a phase factor 
exp(— i2n^). But for a state with exactly one fluxon lo- 
cated next to site i, this factor is exactly cancelled by the 
term exp(i27r^| °r) present in equation (|30|l . leaving 
only the local gauge generator Ui. This shows that a sin- 
gle fluxon moves in this limit as a free quantum particle 
on the dual lattice, with an energy spectrum: 

27T 

e(k) = — cos( — )) — 2A~ 1 (cos(fc E ) + cos(fc„)) 
n 

This spectrum is gapped as long as A remains larger than 
\ c n given by: 

1 



(40) 



u , ,2tt ss 

— = 9(1 -cos — 
AS 4^ v n " 



(41) 
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These equations neglects the renormalization of the 
fluxon spectrum by the fluctuations, but we expect that 
it remains qualitatively correct when these fluctuations 
are taken into account. In particular, even in the regime 
of the strong fluctuations, close to the transition, the 
picture of the fluxons moving with a spectrum e(k) = 
5e + 27^k 2 should remain valid at sufficiently long scales 
as long as A n > A£. For a system with an external charge 
Qi located at site i, we have Ui = exp(i^ L Q i ), so a sin- 
gle fluxon will experience the usual Aharonov-Casher ef- 
fect from this static charge, independently of the value 
of the Chern-Simons coefficient v. Since this interfer- 
ence effect always raises the value of the fluxon energy, 
in comparison to the case Qi = 0, we see that the flux 
attachment mechanism predicted in the continuous £7(1) 
Chern-Simons theory does not operate as long as the sin- 
gle fluxon spectrum remains gapped. In the notations of 
this section (cf. (|IV|) above), this implies that A* = A°. 



VI. CONCLUSION 

We presented (Section II) the general symmetry anal- 
ysis of physical systems with protected degeneracies, i.e. 
with degeneracies that are exponentially weakly affected 
by local perturbations. We have shown that such pro- 
tected degeneracies appear in a system described by a 
wide class of Hamiltonians that commute with two sets 
of integrals of motion, {Pi} and {Qj} but which do not 
commute themselves. These sets of non-commuting op- 
erators should allow for a finite dimensional representa- 
tion. In the simplest case, if PiPj and QkQi commute 
with all other operators the algebra of these operators 
allows two dimensional representation and the states of 
the system are exactly doubly degenerate. For the effect 
of local perturbations to become really small in the ther- 
modynamic limit one also needs that the gap to the low 
energy excitations remains finite. 

We have explicitly constructed a two dimensional lat- 
tice spin model with local interactions and which has 
these integrals of motion. In this model all states are 
exactly doubly degenerate. The behavior of this spin 
model is characterized by the dimensionless parameter, 
J z j J x which physically corresponds to the anisotropy of 
the couplings in different directions. We were able to 
treat it analytically in the regime of large (or small) val- 
ues of this parameter. In this regime the spectrum of the 
system contains 2 L low energy modes where L is the lin- 
ear size of the spin array. The gap between these modes 
and the ground state decreases exponentially with the 
system size. The number of these low energy modes is 
the same as would be the number of edge states but, 
unlike the latter, they are not sensitive to the boundary 
conditions. In order to check the validity of these con- 
clusions for all values of J z /J x we have also performed 
the diagonalization of small arrays (up to 25 spins) and 
concluded that the gap to low energy states remains a de- 
creasing function of the system size for all values but this 



decrease becomes very slow for J z / J x ~ 1. It remains un- 
clear to us, however, whether the system exhibits a new 
phase at these values of the parameter or this apparently 
slower decrease of the gap is a consequence of a critical 
behaviour. 

We have suggested and studied (Section III) two de- 
signs of the Josephson junction arrays and showed that 
their effective low energy Hamiltonians satisfy the sym- 
metry requirements described above and thus their states 
are doubly degenerate and protected from the external 
noise. The simplest of these arrays can be mapped onto a 
spin model with non-local interactions. The non-locality 
of these interactions, however, is not important for the 
protection from the external noise. Further, in these sys- 
tems one can completely eliminate the dangerous low en- 
ergy modes by appropriate boundary conditions. The 
mapping of the Josephson junction array onto a spin sys- 
tem with symmetric Hamiltonian implies that the contin- 
uous superconducting phase can be integrated out. We 
have examined the conditions when this can be done and 
when low energy degrees of freedom corresponding to 
the continuous phase are irrelevant. Summarizing the 
requirements for a physical Josephson junctions we con- 
clude that they are relatively easy to satisfy in medium 
sized arrays (up to 10 x 10 elements) which should be 
quite sufficient to get a noise suppression by 10 orders of 
magnitude. 

The spin models studied in Section II can be also 
mapped onto a discrete Chern-Simons theory on a lattice. 
In order to establish this mapping we have constructed 
(Section IV) a Hamiltonian framework for lattice Chern- 
Simons theories with abelian groups. We argued that, 
in contrast to the continuous theories, such theories gen- 
erally have low energy modes corresponding to the exci- 
tations with large momentum, comparable to the inverse 
lattice spacing. Further, we showed that in a theory with 
a compact group (in particular, in a theory with discrete 
group) the Chern-Simons coefficient, i>, is quantized sim- 
ilar to the quantization of the magnetic flux through the 
torus: v = n/m. In the gauge invariant space of mag- 
netic fluxes the kinetic part of the Hamiltonian of these 
theories can be described as a flux dynamics. Due to the 
presence of Chern-Simons term the motion of fluxes in x 
and z directions does not commute. In the simplest case 
of the Z2, m = 2 theory, the fluxes take only two values 
and their motions in x and z directions anticommute al- 
lowing us to map this theory onto the spin model studied 
in Section II. In a general case the Hamiltonian in flux 
representation becomes very non-local but still this rep- 
resentation is convenient in the limit of large magnetic 
energy when the fluxes are rare. Using this limit, we 
show that the flux attachment to the charge only occurs 
if a single fluxon is gapless. 

The main theoretical issue raised by the Sections IV 
and V of the paper is the precise connection between 
lattice and continuous versions of Chern-Simons theo- 
ries. For the continuous case, there is a sharply iso- 
lated ground-state subspace, whose degeneracy directly 
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reflects the topology of the two-dimensional space on 
which the model is defined. In a recent series of pa- 
per a?^? 1 , various descriptions of these models (in terms 
of wave-functions defined on equivalence classes of loops) 
have been advocated to construct candidate lattice mod- 
els which would exhibit a ground-state sector equivalent 
to a pure (topological) Chern-Simons theory. The ap- 
proach we have followed here starts from a direct quan- 
tization of a lattice Hamiltonian inspired from the con- 
tinuous Chern-Simons theory with an additional kinetic 
term. We found that such construction typically leads 
to the degenerate modes attached to the Brillouin zone 
boundary; it remains to be investigated whether or not 
the presence of these modes spoils the expected proper- 
ties of continuous models (such as statistical transmu- 
tation of external charges). It would be also interest- 
ing to see whether these degeneracies remain for other 
lattice structures. Finally extending this construction 
to non-Abelian discrete groups is clearly desirable, from 
the perspective of enlarging the set of unitary operations 
generated by adiabatic exchanges between charge and/or 
vortex excitations2*32i2i. 
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APPENDIX A: EXACT SOLUTION OF TWO 
CHAIN PROBLEM 

For the diagonalization of the hamiltonian @ it is con- 
venient to rotate the Pauli matrices: t? — > af , rf — > 
erf and use the Jordan- Wigner fermionic representation 
(seeSS) 



(tf = af exp < in ^ a\,a v 




(Al) 

so that © takes the form 

{n-l 
Y,i a t - a 0(«/+i + a i+i) 
i=i 

+2AJ2(4*1- 1/2)1, (A2) 
i=i ) 



where A = J X /2J Z . By means of the Bogoliubov-like 
transformation the hamiltonian (|A2() can be diagonalized 

H = J2E k (2b+b k -l) 

k 

where , bk are the fermionic operators of the eigen- 
modes with eigenenergies 



E(k) = 2J 2 VA 2 + l + 2Acosfc. 

The quasi-continuous spectrum of this hamiltonian can 
be found from the quantization rule for the quasimomen- 
tum k: 



k(n + 1) — arctan 



sin k 



A + cos k 



77171, 



(A3) 



with integer m. For A > 1 this equation has exactly 
n distinct nontrivial solutions, and the set of the corre- 
sponding eigenfunctions is complete. It is not the case 
for A < 1, however. Here the number of the continuous 
spectrum eigenstates is only n — l, so that there should be 
one additional mode - the bound state. To find the latter 
one should look for complex solutions of the dispersion 
equation 1A3(I . Substituting k = tt + vy, we arrive at 



7(n + 1) = - In i - 

n ' 2 I A-e-T 



(A4) 



and, introducing e~ 7 = A + x with small x -C 1, we get 
x m (4- — A) A 2 (™ +1 ), so that and the bound state energy 
is 

2A = 2E{rc + i 7 ) a 4J Z [1 - {J x /2J z ) 2 ]{J x /2J z ) n . 

(A5) 

The formula l|A5jl is valid for (J x /2J z ) n < 1. In the 
fermionic representation the bound-state eigenfunction 
is localized near the ends of the chain within the range 
£ ~ — In A. The localized character of the mode respon- 
sible for the ground state doublet splitting is, however, 
an artifact of the nonlocal representation IjAlfl and, ap- 
parently, does not have much physical meaning. 

A similar result for the splitting can also be obtained 
for an Ising chain with periodic boundary conditions. 
Here the chain does not have ends, and the dispersion 
equation has only solutions, corresponding to the contin- 
uous spectrum. The splitting arises from the following 
effect: the effective boundary conditions for the Jordan- 
Wigner fermions a are periodic or antiperiodic, depend- 
ing on the parity / of their total number (the latter is a 
good quantum number): 



kn 



27T771, 

2tt(to+ 1/2), 



for I 
for I 



1. 



(A6) 



with m — 0, 1, . . . , L — 1. 

The energy of the first excited state is just due to this 
effect; it comes not from any specific single-particle state, 
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but from the entire Fermi sea of the filled energy lev- 
els, each of which is slightly shifted when the parity / is 
changing: 



2?rm\ „ /, 2vr(m + l/2) 
-E \ k= 



(A7) 



Using the Poisson summation formula we arrive at the 
result 



2A * 4.7-4/ 1 (Jx/2Jz) \ j x /2J z r, (A8) 



irn 



which differs from (|A5|I only in the preexponential factor. 



APPENDIX B: THE CRITICAL POINT OF 
THREE CHAIN PROBLEM 

We consider a three-strings ladder with periodic 
boundary conditions along each rung. The corresponding 
Hamiltonian reads 

n— 1 3 n 3 

ff=-^EE^ir^EE^ + i. (si) 

»=1 3=1 *=1 3=1 

where j = 4 is identical to j — 1. We introduce a basis 
of 4 states ip m (with m = 0, 1, 2, 3) on a particular rung 
i, corresponding to the sector with all Pj = 1: 




Then (up to an irrelevant additive constant) the Hamil- 
tonian (IBlll can be rewritten as 



EJ ^J z ^ ^ fimimi+i Jx ^ ^ r ? ;, 



where the matrix T 



r = 



1 o 
i i 

V i i 



1 1 \ 
i i 

o 1 



plays a role, similar to that of the a^-operator for the 
two chain problem. This is the hamiltonian of the one 
dimensional q = 4 Potts model in a " transverse field" . 

Consider now an asymmetric classic two-dimensional 
q = 4 Potts model with the Hamiltonian 

f3H~2d = ~K Z ^ S m ik m i+lk -K x ^2 <5m sfc m, fc + 1 , (B2) 



ik 



ik 



The transfer-matrix for this system (in the fc-direction) 



is 



{n-l 
K z <W +1 j> Y[ (le K * + f i) , (B3) 



i=l 



where 1 is 4 x 4 unity matrix. Using the identity 



r =i{(e 3 ' l + 3e-' I ) + ( e 3 ' I - e -' l )f} 



for matrix T, we can rewrite (IB3I) in the form 



f n— 1 n ^ 

f = C exp i K z ^ 6 



(B4) 



where C is an irrelevant constant and h is determined by 
e K * = (e 4h + 3)/{e 4h ~ 1). 

The line of critical points for the asymmetric two- 
dimensional q = 4 Potts model is governed by the relation 
(seeii) 

(e*"-l)(e*--l) = 4; 

in terms of _fT z , /i this relation takes simple form K z = Ah. 
On the other hand, the matrix l|B4|) describes the time 
evolution of the quantum system with the Hamiltonian 
(IB!)) and with 4J Z = , J x = -j^, where At — > is 
an infinitesimal time-interval. Thus, we conclude, that 
the quantum phase transition in our initial three chain 
system takes place at the symmetric point J x = J z . 

Unfortunately, the solution of the two-dimensional q = 
4 Potts model away from the critical line is not known 
and, in contrast to exactly solvable two-chain model, we 
can not find the dependence of the gap A on the param- 
eter J x / J z in the full range of this parameter. 
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